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ABSTRACT 



o: 

(N. Context. It is believed tiiat the majority of stars form in clusters. Therefore it is likely that the gas physical conditions 

1— ( ■ that prevail in forming clusters largely determine the properties of stars that form in particular, the initial mass function. 

^ I Aims. We develop an analytical model to account for the formation of low-mass clusters and the formation of stars 

' ' , within clusters. 

^\ • Methods. The formation of clusters is determined by an accretion rate, the virial equilibrium as well as energy and 

^SJ ' thermal balance. For this latter, both molecular and dust cooling are considered using published rates. The star 

I distribution is computed within the cluster using the physical conditions inferred from this model and the Hennebelle 

I I, & Chabrier theory. 

■ Results. Our model reproduces well the mass-size relation of low mass clusters (up to few ~ 10^ Mq of stars corresponding 

Oto about five times more gas) and an initial mass function that is i) very close to the Chabrier IMF, ii) weakly dependent 
on the mass of the clusters, in) relatively robust to (i.e. not too steeply dependent on) variations of physical quantities 

(~| , such as accretion rate, radiation, and cosmic ray abundances. 

^ ^ Conclusions. The weak dependence of the mass distribution of stars on the cluster mass results from the compensation 

I ■ between varying clusters densities, velocity dispersions, and temperatures that are all inferred from first physical 

O ' principles. This constitutes a possible explanation for the apparent universality of the IMF within the Galaxy although 

, variations with the local conditions may certainly be observed. 
C/5 . 

, Key words. Instabilities - Interstellar medium: kinematics and dynamics - structure - clouds - Star: formation - 

I Galaxies: clusters 

T— I ■ 

■ 1. Introduction proportional to p^^^. While this later idea is interesting, it 

J ' raises a few questions. First, the first generation of stars is, 

J — . Since the pioneering vi^ork of Salpeter (1955), the origin of at least at the beginning of the process, not influenced by 

\^ ' the initial mass function (IMF; e.g. Kroupa 2002, Chabrier the radiative feedback. Second, even when the stars start 

2003) remains one of the most fundamental questions in forming, the regions in which heating is important remain 

astrophysics. More particularly, its apparent universality, limited to the neighborhood of the protostars. Thus it is 

■ Bastian et al. (2010), remains debated and challenging, not yet demonstrated that most stars will be affected by 
Indeed, as described by these authors, the IMF has now this effect and that it is sufficient to make the IMF uni- 
been determined in many environments, and although some versal. Moreover the recent simulations by Krumholz et al. 
of the results may be interpreted as evidence of variations, (2012) suggest that indeed the gas temperature is much 
neither systematic, nor undisputable evidence for such vari- too high when radiative feedback is included, leading to 
ations of the IMF are unambiguously reported. Therefore an IMF that is inconsistently shifted toward high masses. 

] even if, as it will probably eventually turn out, some varia- On the contrary when outflows are included, the radiation 

tions are finally clearly established, in many environments, can escape along the outfiow cavities and they obtain IMF 

the variations of the IMF remains limited. which are close to the observed one. 

Various theories have been proposed to explain the ap- Another type of arguments invokes the variation of the 

parent constancy of the IMF. These theories often rely on effective polytropic index, 7, which in particular presents a 

the independence of the Jeans mass on the density. For local minimum at about 10^ cm~'^ (e.g. Larson 1985) due 

example, Elmegreen et al. (2008) computed the gas tem- to the transition between cooling dominated lines and dust, 

perature in various environments and found a very weak This in particular has been proposed by Bonnell et al. (2006) 

dependence of the Jeans mass on the density since the lat- and Jappsen et al. (2005). However, it remains unclear that 

ter increases with temperature. Bate (2009) and Krumholz it is actually the case because the various simulations did 

(2011) estimated the temperature in a massive collapsing not clearly establish that changing the cloud initial condi- 

clump in which the gas is heated by the radiative feedback tions while keeping a fixed equation of state would lead to 

due to accretion onto the protostars. Bate (2009) found that a CMF peaking at the same mass. Moreover Hennebelle & 

in the vicinity of the protostars, the Jeans mass is typically Chabrier (2009) compare their predictions with the peak 
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position found in numerical simulations of Jappscn ct al. 
(2005) and find a good agreement. Yet in the Hennebelle 
& Chabrier theory, a local minimum of 7 does not deter- 
mine the peak of the CMF which still depends on the Mach 
number for example. 

More generally, these theoretical arguments assume that 
the Jeans mass is the only parameter that determines the 
IMF. This sounds rather unlikely because a distribution like 
the IMF is not entirely determined by a single parameter 
(peak position, width, and slope at high masses) , moreover, 
analytical theories like the one proposed by Hennebelle & 
Chabrier (2008, HC2008) and Hopkins (2012) explicitly de- 
pend on the Mach number. Although no systematic explo- 
ration of the Mach number influence on the core; mass func- 
tion has been performed in numerical simulations, Schmidt 
et al. (2010) have explored the role of the forcing of the 
turbulence. They showed in particular that the core mass 
function (CMF) is quite different when the forcing is ap- 
plied in pure compressible or pure solcnoidal modes. This 
clearly suggests that the CMF is affected by the veloc- 
ity field. Even more quantitatively, Schmidt et al. (2010) 
found a good agreement between the CMF they measure 
and the analytical model of HC2008 which seemingly con- 
firms the Mach number dependence of this model. From a 
physical point of view, it is well established that the density 
probability distribution function (PDF) is strongly related 
to the Mach number (Vazqucz-Scmadcni 1994, Padoan et 
al. 1997, Passot & Vazquez-Semadeni 1998, Kim & Ryu 
2005, Kritsuk ct al. 2007, Federrath et al. 2008, Audit' & 
Hennebelle 2010). Accordingly, because the density PDF is 
clearly important with regard to the Jeans mass distribu- 
tion within the cloud, it would be quite surprising if the 
Mach number had no influence on the CMF. 

Another line of explanation regarding the universality 
of the IMF has been proposed by HC2008, who argue that 
because the peak position of the IMF is proportional to the 
mean Jeans mass and inversely proportional to the Mach 
number, there is a compensation because from Larson re- 
lations, the density decreases when the velocity dispersion 
increases and thus the mean Jeans mass and the Mach num- 
ber increase at the same time (see Eq. 47 of HC2008 and 
figure 8 of HC2009). One problem of this explanation is, 
however, that Larson relations present a high dispersion 
and that there are clouds for example, with the same den- 
sity but different Mach numbers that may therefore lead to 
different IMF in particular at low masses. 

Generally speaking, all approaches that attempt to un- 
derstand the universality of the IMF suffer from the vari- 
ability of the star- forming cloud conditions. This clearly 
emphasizes the need for a better understanding of the phys- 
ical conditions under which stars form. In this respect, an 
important key is that most stars (say 50-70%) seem to form 
in clusters (Lada & Lada 2003, Allen et al. 2007), which is a 
strong motivation to study the formation of clusters. Note, 
however, that Bressert et al. (2010) moderate this picture 
to some extent. 

A word of caution is nonetheless necessary here. The 
constancy of the IMF is debated and some observations 
even in the Galaxy may indicate that some variability 
has already been observed (see e.g. Cappellari et al. 2011 
for early types galaxies). For a discussion and good sum- 
mary on this problem we refer the reader to Dib et al. 
(2010). Moreover, these authors have developed an analyt- 
ical model of the core formation within proto-clusters that 



takes into account the turbulent formation of dense cores 
as well as their accretion of gas that could modify the CMF 
and lead to IMF variability. Along a similar line, Dib et al. 
(2007) showed how the coalescence of cores can explain the 
development of a top-heavy IMF. 

In this paper, we first develop an analytical model for 
the formation of clusters, more precisely, the formation of 
proto-clusters, i.e. the gas dominated phase that eventually 
leads to star-dominated clusters. Our model relies on the 
gas accretion onto the proto-clusters from parent clumps 
that feed them in mass and in energy. This allows us to 
predict the physical quantities such as radius, mean den- 
sity, and velocity dispersion within the proto-clusters as a 
function of the accretion rate. By comparing with the data 
of embedded clusters from Lada & Lada (2003), we can esti- 
mate the accretion rate onto these proto-clusters and verify 
that our model fits the observational data well. In a second 
step, we calculate the gas temperature within the proto- 
clusters by computing the various heating and cooling, and 
we apply a time-dependent version of the model of HC2008 
to predict the mass spectrum of the self-gravitating con- 
densations and in particular to study their variability with 
proto-cluster masses and accretion rates. 

The second part of the paper presents the analytical 
model of the proto-clusters and the comparison with the 
observational data. The third part is devoted to the calcu- 
lation of the thermal balance as well as to the description 
of the HC2008 theory. In the fourth part, we calculate the 
mass spectra within the proto-clusters. The fifth section 
concludes the paper. 



2. Analytical model for low mass cluster formation 

Clusters are likely to play an important role for the forma- 
tion of stars in the Milky Way (e.g. Lada & Lada 2003) 
and probably for most galaxies. It seems therefore a ne- 
cessity, in order to understand how star forms to obtain 
a good description of the gas physical conditions within 
proto-clusters. For that purpose, we develop here an ana- 
lytical model that is based on the following general ideas. 
First, proto-clusters are initially likely to be gravitationally 
bound entities. However, proto-clusters are not collapsing, 
which means that a support is actually compensating grav- 
ity, which we assume; is the turbulent dispersion. This will 
be expressed by applying the virial theorem to the proto- 
cluster. Because turbulence is continuously decaying with 
a characteristic time on the order of the crossing time, it 
must be continuously sustained. We assume that the con- 
tinuous accretion of gas into the proto-cluster from the par- 
ent clump is the source of energy. Indeed, accretion-driven 
turbulence has been recently proposed to be at play in var- 
ious contexts (Klessen & Hennebelle 2010, Goldbaum et 
al. 2011). Our first step is therefore to discuss the par- 
ent clumps and the resulting accretion rate onto the proto- 
cluster. Once an accretion rate is inferred, an energy bal- 
ance can be written and together with the relation obtained 
from the virial theorem lead to a link between the mass, the 
radius, and the velocity dispersion within the cluster. 
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2.1. Parent cloud, protocluster and accretion rate 
2.1.1. Definitions and assumptions 

Let us consider a clump of mass Mc ~ 10'^ — IO^Mq 
and radius Rc that follows Larson relations (Larson 1981, 
Falgarone et al. 2004, Falgarone et al. 2009): 



no 



Rc 
Ipc 



Rc ^ 



(1) 



where Uc is the clump gas density and arms the internal 
rms velocity. The exact values of the various coefficients 
remain somewhat uncertain. Originally, Larson (1981) es- 
timate r]d — 1.1 and rj ~ 0.38, but more recent estimates 
(Falgarone et al. 2004, 2009) using larger sets of data sug- 
gest that r]d ~ 0.7 and rj ~ 0.5. These later values agree weh 
with the estimate from numerical simulations of supersonic 
turbulence (e.g. Kritsuk et al. 2007) and we therefore use 
them throughout the paper although we will compare the 
accretion rates obtained for both sets of values. 

As explained above, we consider the formation of a clus- 
ter within the parent clump. Because we are essentially 
considering the early phase during which the gas still dom- 
inates over the stars, we refer to it as the proto-cluster. Let 
M* be its mass and a^.^c = Mc/M^, be the ratio between the 
parent clump and proto-cluster masses. Let us stress that 
M» is the mass of gas within the proto-clusters and not the 
mass of stars. Throughout this work the star component is 
not considered. 

Obviously, the proto-cluster is accreting from the par- 
ent clump, which is therefore regulating the rate at which 
matter is delivered onto the proto-cluster. The immediate 
question concerns then the value of the accretion rate ? To 
answer this we consider two different approaches. First we 
assume that the proto-cluster clump is undergoing Bondi- 
type accretion, i.e. accretion regulated by its own gravity. 
However because this process can be sustained only if suf- 
ficient matter is available for accretion, we also explore the 
possibility that the accretion onto the proto-cluster is regu- 
lated by the accretion onto the parent clumps. In this case, 
the accretion is due to the mechanism that is at the ori- 
gin of the clump formation and that also sets the Larson 
relations. Both assumptions lead to similar, although not 
identical values and dependence. 

2.1.2. Bondi-type accretion rate 

The Bondi-Hoyle accretion rate that an object of mass 
is experiencing in a medium of density p and sound speed Cs 
is given by M = An pG'^ / (Bondi 1941). In a turbulent 
medium the accretion rate is not well established. While a 
simple expression is given by M = in pG'^ / {c^ + a'^)^/^, 
where a is the velocity dispersion, Krumholz et al. (2005) 
obtain and discuss more refined estimates. We stress that 
these estimates assume that the accreting reservoir is infi- 
nite, which, as already mentioned may not be true in the 
present context. Moreover, these accretion rates are valid 
only for point masses and it is unclear whether they apply 
to more extended objects such as proto-clusters. 

To estimate an accretion rate, here we used the modified 
Bondi expression. The density p and the velocity disper- 
sion a are those from the parent clump stated in Eqs. ([T]). 
Because these quantities depend on the parent clump mass, 
Mc, and because likely increases with Afc, it is clear that 
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Fig. 1. Accretion rate onto the proto-cluster against clump 
mass. BLl and BL2 represent Bondi-Lyttleton accretion 
for T]d = 1.1, r\ — 0.4 and r\d = 0.7, r\ = 0.45, respectively. 
Tl and T2 represent the turbulent accretion rate with the 
same values. Upper (lower) BL curves are for parent clump 
over proto-cluster mass of 3 (5). 



the dependence of the accretion rate is not proportional to 
Ml but rather has a more shallow dependence essentially 
because p and a decreases and increases respectively when 
Mc increases. To arrive at a more quantitative estimate, we 
assume that a,.c does not depend on Mc and is therefore 
constant. We obtain 



M = 



inpG'^Ml 



AnMpTioG'^Ml 



a2ca3(i?/lpc)')^+3') 



cx a. 



3(2~(T;rf+r,)) 

-2^ is-va) 



(2) 



For rid = 0.7, rj = 0.45, we get M cx M while for 
rjd = 1.1, rj — 0.4, we obtain M cx Af^*'-^. Figure [1] shows 
the resulting accretion rate, BLl and BL2 referring to the 
first and second sets of parameters respectively, the values 
a*,c = 3 (upper curves) and a*_c = 5 (lower curves) were 
used. This likely corresponds to upper values of the accre- 
tion rate because higher values of a*,c lead to lower values 
of M, while for lower values of a*_c, the mass within the 
parent clumps is equal to or even smaller than, the mass 
within the protocluster and it is therefore extremely un- 
likely that Bondi-type accretion is relevant at all. 

2.1.3. Turbulent accretion rate 

The second type of accretion we consider is the one that 
is responsible for the clump formation. The exact nature 
of these clumps is still a matter of debate but simulations 
have been reasonably successful to explain at least the low- 
mass part of the distribution (e.g. Hennebelle & Audit 2007, 
Banerjee et al. 2009, Klessen & Hennebelle 2010) in which 
they can essentially be seen as turbulent fluctuations. In 
this view, as the properties of the clumps stated by Eqs. ([T]) 
are the result of the turbulent cascade, that is to say, the 
material within the clumps is delivered by the compressive 
motions of the surrounding diffuse gas, it seems reasonable 
to construct an accretion rate out of the Larson relations. 
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The clump crossing time, Tc is about 2R/ [a / ^/i) where 
the factor ^/i accounts for the fact that it is the one- 
dimensional velocity dispersion that is relevant for estimat- 
ing the crossing time. 

The accretion rate of diffuse gas onto the clump is ex- 
pected to be 



M 



Mr 



Mr 



Tc 2i?,/((T™,/V3)' 



°- ? 



Mc ^""^ ^ (^y nomp ) (1 pc) , 

1714 I 



(3) 



9 X IQ-'^Mq yr~i 



f Mr \ ^-"d 



V IO^Mq 
(To / rifi 



0.8kms-i VlOOOcm-3 



1-'? 



where mp is the mass per particle. Note that to estimate 
7714 we used rjd — 0.7 and rj — 0.45. 

The Larson relations that describe the mean properties 
of the clump in which the cluster forms are likely a direct 
consequence of the interstellar turbulence. Thus, the accre- 
tion is likely to last a few clump crossing times, which is 
typically the correlation time in turbulence. 

If the clump is sufficiently dense, it undergoes gravita- 
tional collapse and forms a cluster of mass . Here we at- 
tempt to understand the cluster formation phase. For that 
purpose, we assume that: 

- as already explained it starts with a gas-dominated phase, 

i.e. the mass of the cluster is dominated by the mass of 
its gas and not its stars 

- this phase is quasi-stationary, i.e. the cluster properties 

evolve slowly with respect to its crossing time. 

- the accretion rate onto the cluster is comparable to the 

accretion rate of diffuse gas onto the clump and is there- 
fore close to the value given by Eq. 

Of these three assumptions, the second one is certainly the 
least obvious. Indeed, initially no star will have formed by 
definition, while what we assume in practice regarding the 
accretion rate is simply that it scales as M oc AI^"^" and 
we fit the coefhcient using observed proto-clusters. Taking 
time-dependence into account would certainly be highly de- 
sirable at some stage but will probably not modify the re- 
sults very significantly. 

It is useful to express the accretion rate as a function of 
the proto-cluster mass and we write 



— TTlA I - 



M — 7124 



(4) 



Assuming that the parameter rh^ stays constant, this im- 
plies that the accretion rate onto the proto-clusters is con- 
trolled by Mc, the mass of their parent clumps, and thus by 
a*,c = Mc/M*, which effectively quantifies the strength of 
the accretion onto the proto-cluster. In practice, 7714 could 
vary while q;*^c would be constant, it is equivalent, how- 
ever, and does not make any difference if one uses one or 
the other. Below we use a»^c to quantify the accretion onto 
the proto-cluster. 



2.2. Virial equilibrium 

Because a cluster initially is a bound system, which unlike 
a protostellar core is not globally collapsing, it seems clear 
that on large scale some sort of mechanical equilibrium is 
established, to describe which we use the virial theorem. 
When applying the virial theorem to the cluster, it must 
be taken into account that it is gaining mass by accretion 
and is therefore not a close system. In Appendix A, we show 
that when applied to an accreting system, the expression 
of the virial theorem becomes 



dMl 



0. 



(5) 



In this expression, is the cluster mass, cr* is the rms 
velocity of the gas within the cluster, = 47r/3i?2 is its 
volume, and Pram is the external pressure exerted by the in- 
falling clump gas onto the cluster which dominates over the 
thermal pressure. Assuming that the gas within the parent 
clumps is gravitionally attracted by the proto-cluster, the 
infall velocity is simply the gravitational freefall, and we 
obtain 



2GM^ 



(6) 



The ram pressure is given by Pram = Pinfi'hif while the 
accretion rate M leads to the relation M = 47rP^pinfWinf , 
therefore leading to 



p — 

ram — 



M 
47rP2 



Vinf 



2V2ttR 



5/2 



(7) 



The gravitational energy of a uniform density sphere is 
well-known to be -{3/5)GM^/R^, and thus we obtain 



,V2 dR: 

M^ 



dM^ 



M,ai - SPramK 



3GMl 
5 i?. 



0. 



(8) 



2.3. Energy balance 

As stated in Eq. ([5]) , the proto-cluster is confined by its own 
gravity and by the ram pressure of the incoming flow. The 
velocity dispersion of the gas it contains resists these two 
confining agents. To estimate its magnitude, we assume an 
equilibrium between energy injection and turbulent dissi- 
pation 



2Tc.r.t 



(9) 



where Tcct is the cluster crossing time and is about Tcct — 
2R„/a^,/ Vs. -Bint and E'cxt are the external and internal 
source of energy injection which compensate the energy 
dissipation stated by the left side. Note that in principle 
a complete energy equation could be written that would 
entail thermal energy but also terms taking into account 
the expansion or contraction of the proto-cluster (see e.g. 
Goldbaum et al. 2011). However, the exact amount of en- 
ergy dissipated by turbulence is not known to better than a 
factor of a few and the same is true regarding the efficiency 
with which the accreting gas can sustain turbulence (giv- 
ing that a certain fraction can quickly dissipate in shocks) . 
Since these two contributions are dominant, we consider a 
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simple balance at this stage from "which physical insight can 
be gained. 

From Eq. we obtain 



dM/dt ocM" ' 



2/3 



int 



(10) 



Note that in this work we do not consider any internal 
energy source such as supernova explosions, jets, winds, and 
ionizing radiation from massive stars, accordingly Ei^t ~ 
is assumed in this paper. 

The external source of energy is generated by the accre- 
tion process itself because the kinetic energy of the infalling 
material triggers motion within the cluster. This energy flux 
is on the order of {l/2)Mvf^f ~ GMM^/R^. It is slightly 
larger than this value, however, because this expression cor- 
responds to the kinetic energy of the gas when it reaches 
the cluster boundary. In practice the gas accreted onto the 
cluster continues to fall inside the cluster and gain addi- 
tional energy. In appendix B, we infer the corresponding 
value, whose expression is 



E,- 



6 GMM, 
5' 



GMlVh 



(11) 



Thus with Eq. ©, ^ and ([U]) we get 



3 GM, 
5 i?. 
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1 dRl 



(12) 



where K = 12\/3/5 although as discussed above it suffers 
from large uncertainties. 

2.4. Result: Mass vs radius of clusters 
2.4.1. Simplified expression 

Before deriving exact solutions of Eq. (fT2)) and to obtain 
some physical hint, we start by discussing the simplified 
case where the terms proportional to Fcxt and are ig- 
nored. In this case we simply have 



4.3 



(gm) 



2/3 



(13) 



Using Eq. ([3]), we obtain 



4.3T^ 



G-i(m4)'^ " 



V lO^Mp 



2riacc\ 3-2170 



X (i?^)l/(l-2..ee/3)^ 



(14) 



where we recall that q;,,c — AIc/AI^, is the ratio of cluster 
over clump masses. As seen from Eq. ([3]), rjacc — 0.75 (rj — 
0.4 — 0.5), therefore 1/(1 — 27yacc/3) — 2, which is the value 
that we adopt to perform this simplified calculation. We 
obtain 



~ 4.32G"2/3(^4)4/3(^^ JlO4M0)i?2, 

which implies 

4/3 

M* ~ 730 Mq X a*,c 
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Fig. 2. Upper panel: Mass of the cluster as a function of 
its radius as predicted by Eq. ([T^ . i.e. for a turbulent type 
accretion for which M oc M° ''^. Three values of a^,^c are 
displayed. The stars correspond to the embedded clusters 
listed in table 1 of Lada & Lada (2003) for which the mass 
has been multiplied by 5 to account for a star formation 
efficiency of about 20%. Lower panel: same as lower panel 
for M c>c Af . The thick solid line corresponds to the solid 
line of the upper panel for comparison. 



As suggested by Eq. ([3]), the value of m4 is expected 
to be about a few lO~'^M0yr~^ while assuming that the 
mass of proto-cluster is comparable but smaller than the 
mass of its parent clump, a,^c which is likely to be about 
1-4. Accordingly, for a 1 pc cluster, one typically expects a 
mass of about 700-2000 Mq of gas. 



(15) 2.4.2. Complete expression for turbulent-type accretion 



We now solve Eq. (jl2l) in the general case, i.e. without 
neglecting some of the terms. It can be shown that assuming 
Vacc = 0.75, the solution of this equation can still be written 
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as M* = Af°(i?,/lpc)2., where M° satisfies 



2.4.4. Comparison with observations of embedded clusters 
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(17) 



Solving this equation numerically for a^.c^l, 2 and 4, we 
infer M° ~ 4144, 1381 and 414 M©, respectively, which is 
about 1.8 times lower than the value estimated in Eq. ((T6)) . 

We note that with oc we derive that the gas den- 
sity and the column density follow cx R~^, respectively, 
while E is constant. Quantitatively, we have the three rela- 
tions 



S,. = 2.8 X 10-cm-^, 



(18) 



where Af° is defined by Eq. = A/*/(47r/3mpi?3) 

and S,^c = 2n^,R^,. 

Finally, from Eqs. ([3]), ([TU| and (ITT|) . we obtain the in- 
ternal gas velocity dispersion within the proto-clusters 



Vlpc 



1/2 



(2V3G)- 



1/3^1/3 ^a^^y^Vi?* ^ 

104 Af©; ^ipc^ 



2.4.3. Bondi-type accretion 
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As discussed above, combining Bondi-type accretion and 
Larson's relations, we infer accretion rates whose mass de- 
pendence follows dM/dt oc M^'^"" with rjacc — 1- Using 
this canonical value, it is easy to see that the solution of 
Eq. ini) is = Af,"(i?/(lpc))3, where M° satisfies the 
following equation: 



^GAf," + V2Gi/2g;(A/0)V2 = 

2/3 



1' • I 



\ 104A/(r 



0.75 



X (lpc)3/2. 



For canonical values of a*_c = 2 and 7714 = 10 M© s ^, 
the relation M» ~ 1380 M© (i?*/lpc)2 holds. To test our 
model, we compare it with the data of embedded clusters 
reported in table 1 of Lada & Lada (2003). Because the 
mass quoted in this table is the mass of the stars, one must 
apply a correction factor to obtain the mass of the gas. 
In table 2, Lada &: Lada give an estimate of this ratio for 
seven clusters. The star mass over gas mass ratio is typically 
between 0.1 and 0.3 with an average value that is about 0.2. 
Thus to perform our comparison, we simply multiply the 
star masses of table 1 by a factor of 5. By considering a 
unique star formation efficiency, we certainly increase the 
dispersion in the data. However since a few values for the 
star formation efficiency are available, this is unavoidable. 

The upper panel of Fig.[2]shows the results for three val- 
ues of Q!*_cj the clump over proto-cluster mass ratio, namely 
1, 2 and 4. Evidently, a good agreement is obtained with 
q;*.c — 2 over almost 2 orders of magnitude in mass, al- 
though the dispersion of the observed values is not negli- 
gible (factor 2). This may suggest that the relation stated 
by Eq. ^ is relatively uniform throughout the Milky Way. 
Because this is likely a consequence of interstellar turbu- 
lence, it may reflect the universality of its properties. Note 
that as emphasized in Murray (2009), high- mass clusters 
(> 10^ Mq) present a different mass-size relation. This may 
indicate that for more massive clusters other energy sources 
than stellar feedback should be considered. It is also likely 
that the feedback strongly affects the proto-clusters when 
enough stars have formed. In particular, gas expulsion is 
likely to occur and feedback probably sets at least in part, 
the star formation efficiency (e.g. Dib et al. 2011). This late 
phase of evolution is not addressed here, however. 

The lower panel of Fig. [5] shows the results for the 
Bondi-type accretion, i.e. Af oc A^ as emphasized in sec- 
tion [^321 The thick solid line reproduces the solid line of 
the upper panel to facilitate comparison. Clearly the agree- 
ment is not as good as for the turbulent-type accretion for 
which M oc Af°-^^, which is clearly in favor of this latter 
case. Therefore we restrict ourselves to this latter case in 
the following. Physically, it suggests that indeed large scale 
turbulent fiuctuations may be regulating the accretion onto 
the proto-stellar clusters. 

To summarize, the good agreement we obtain between 
our theory of proto-cluster formation and the data of em- 
bedded clusters suggests that i) accretion-driven turbulence 
(Klessen & Hennebelle 2010, Goldbaum et al. 2011) is at 
play in low-mass proto-clusters, ii) the accretion rate onto 
proto-clusters is reasonably described by Eq. Q with a 
value of a*^c — 2 and M oc Af° ''^. Therefore we adopt 
these values as fiducial parameters for the remainder of the 
paper. 



(20) 3. Thermal balance and mass distribution 



(19) 



That is to say, a dependence of the accretion rate M oc M 
leads for the protoclusters to a mass-size relation M oc R^, 
i.e. density is independent of the mass. 



In this section we compute the temperature of the gas 
within the cluster and we also recall the principle and the 
expression of the Hennebelle & Chabrier theory considering 
a general equation of state. 
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3.1. Heating rate 



3.2. Cooling rate 



3.1.1. Heating by turbulent dissipation 



As discussed in section [^31 the turbulent energy within the 
cluster is continuously maintained by the accretion energy. 
The turbulent energy eventually dissipates, being converted 
into thermal energy, which is then radiated away. The gas 
within the cluster is thus subject to a mechanical heating 
equal to the expression stated by Eq. ([TT|) . The mechanical 
heating per particle in the cluster is given by 



_ G GMrrip 

^ turb — "Z" T~t 

5 it. 



ergs 



(21) 



Note that as for the turbulent energy balance, a complete 
heat equation could be written, but as discussed above, 
large uncertainties hampered energy dissipation. Moreover, 
this heating represents an average quantity but may greatly 
vary through space and time in particular because turbu- 
lence is intermittent. Indeed, one may wonder whether the 
amount of energy dissipated per unit of time could not de- 
pend on the gas density. Various authors (Kritsuk et al. 
2007, Federrath et al. 2010) found a weak dependence of 
the Mach number on the density, the former decreasing as 
the latter increases. However, the dependence is extremely 
weak, Ai cx p~^ '^^. It seems therefore a reasonable as- 
sumption to treat the mechanical heating as being uniform 
throughout the cluster. 

With the help of the results of the preceding section, 
we can estimate this heating. For a cluster of radius ~ 1 
pc, the mass is about 10'^ M© and the accretion rate 
(2 X lOVlO'')" '^'^ X 10-3 Mq yr-i where a^.c = 2 has been 
assumed. Thus Fturb — 1-6 x 10^^'' erg . 

3.1.2. Cosmic ray heating 

In this work, wc assume that the proto-cluster is embedded 
in its parent clump, whose mass is a few times larger. Since 
the total gas mass of the proto-cluster is 100 to 10* Mq, 
the mass of the parent clump is typically a few times this 
value, as discussed above. Hence the column density of the 
parent clump is about 



a. 



1/2.3 



(4V3)nomp(lpc)o^ 

2i?cnc = 2no(lpc)°-^i?°-3 ~ 5 - 10 X lO^^ cm^^ (22) 

This implies that the visual extinction of the gas surround- 
ing the proto-cluster is typically a few Ay. In the same 
way, it is typically about 10 or more throughout the clus- 
ter, as stated by Eq. ([T5)) . This means that it is fair to 
consider that proto-clusters are sufficiently embedded and 
optically thick to neglect the external UV heating. On the 
other-hand, cosmic rays are able to penetrate even into well 
shielded clouds, providing a heating rate of 



10 



-27 



c 



3 X lO-i^s 



ergs 



(23) 



where ^ is the mean cosmic ray ionization (e.g. Goldsmith 
2001). Note that our notation here differs from the choice 
that is often made because F is the heating per gas particle 
and not a volumetric heating. 

Comparison between this value and Fturb reveals that 
they are indeed comparable, with the higher dominating 
the former in massive proto-clusters. 



In the proto-cluster, the mean density is typically a few 
1000 cm^^ and the visual extinction is about 10. In these 
conditions, the gas is entirely molecular and well screened 
from the UV, as already discussed. Two types of cooling 
processes must be considered, the molecular line cooling 
and the cooling by dust. 

3.2.1. Molecular cooling 

For the molecular cooling we use the tabulated values cal- 
culated and kindly provided by Neufeld et al. (1995). These 
calculations assume that the gas is entirely screened from 
the UV background, which is the case here. It is also as- 
sumed that the linewidth is the same for all species and 
is dominated by microturbulence. The corresponding val- 
ues of the resulting cooling are displayed in Figure 3a-3d of 
Neufeld et al. (1995). The table provided has temperatures 
between 10 and 3000 K, densities between 1000 and 10^° 
cm"'^ with logarithmic increments. The column density per 
km s^^ is equal to 10^°, 10^^ or 10^^ cm"^. The cooling per 
particles can typically be written as 



A. 



mol 



Aon^T''ergs-\ 



(24) 



where c is small and d about 2 — 2.5 (see also Goldsmith 
2001 and Juvela et al. 2001). Note that A/n represents the 
cooling per particle. 

3.2.2. Dust cooling 

The dust must also be taken into account in the thermal 
balance of the gas. The amount of energy exchanged per 
unit of time between a gas particle and the dust is (e.g. 
Burke & Hollenbach 1983, Goldsmith 2001) 



A 



dust 



2 X IQ-^-'niT - T, 



d) 



T 

n ■ V lOK 

where Td is the dust temperature. 

3.2.3. Dust temperature 



ergs 



(25) 



As shown by Eq. (|25|) . it is necessary to know the dust 
temperature to compute Adust- For this purpose we closely 
follow the work of Zucconi et al. (2001). The dust temper- 
ature is the result of a balance between the dust emission 
and the absorption by the dust of the external infrared ra- 
diation, leading to 



QyBy{Td{r))dv = i QyJ^{r)dv, 
Jo 



(26) 



where is the grain absorption coefficient, is the 
Planck function and </i/(r) is the incident radiation field 
given by 



Mr) 



JIS 

An 



exp(-T^(r, 6^, 0))drj, 



(27) 



where J*'* is the interstellar radiation and Tj, the optical 
depth. Equation (|27)) represents the interstellar radiation 
attenuated by the dust distribution within the proto-cluster 
and the parent clump. 
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The values of Qu and J^f are given in Appendix B 
of Zucconi et al. (2001) and are used here. To calcu- 
late Jv{t) it is necessary to specify the spatial distribu- 
tion of the dust. We assume that the proto-cluster is em- 
bedded into a spherical clump of mass Mc = Q!*_c-^^*, 
the column density Sc = 2Rcnc, can be estimated as in 
Eq. (j22p . We make the simplifying assumption that the ex- 
ternal radiation that reaches the edge of the proto-cluster 
is J,y(0) = J** exp(— Sc/2 X Qy). Hence because the cluster 
is assumed to be on average uniform in density, the radia- 
tion field within the cluster can be estimated as described 
in Appendix A of Zucconi et al. (2001), that is 



Mr) 



MO) 
47r 



exp(-T^(r, 0))27rsin(6l)d6l, 



TAr,d) - Q.S:(^l-(r/i?,)2sin2((?)- (r/i?,) cos(0))(28) 

where is the proto-cluster column density toward the 
center. Equation ([^5]) represents the integration of the ra- 
diative transfer equation in all directions through a cloud 
of radius i?* that has a uniform density. 

To find the dust temperature at a given position within 
the proto-cluster, we solve Eq. (|26l) using the wavelength 
dependence of J** and Qi, given in Appendix B of Zucconi 
et al. These values represent a fit from the interstellar 
radiation field given by Black (1994) and grain opacities 
from Ossenkopf & Henning (1994), respectively. For the 
latter, standard grain abundances are assumed. In prac- 
tice this entails the calculations of integrals of the type 
indicated in Eq. (9) of Zucconi et al. Because we do not 
only calculate the dust temperature at the proto-cluster 
centre, Eq. (10) of Zucconi et al. must be replaced by 
J^2l' ty{e^l>{t/l3,) - l)Jt{r)dt. Since Jt needs itself the 
calculation of an integral (that we express as Eq. (A. 3) of 
Zucconi et al.) as shown by Eq. ([28|. the calculations re- 
quire some integration time. We verify that our results com- 
pare well to the results shown in Fig. 2 of Zucconi et al. and 
we obtain temperatures of about ~ 8 K in the proto-cluster 
center and 9.5 K at the edge. 

Finally, because we ought to determine a single dust 
temperature within the proto-cluster, we compute the mean 

dust temperature as fd = /^f * Td{r)'iTTr'^dr /{Att /ZR 
which is typically equal to about 9 K. 



^) 



3.3. Temperature distributions 

To find the temperature of a proto-cluster whose density 
and column density per km s"""^ are known, we must solve 
for the thermal equilibrium 



turb 



=) = A 



mol 



■A 



dust ■ 



(29) 



We start by computing Tdust as described above. To 
solve Eq. (j29p . we use an iterative method, performing a 
first order interpolation in logarithmic density, column den- 
sity per km s^^ and temperature to derive the tabulated 
molecular cooling. Note in passing that the column den- 
sity per km s~^, which is equal to E*/(cr*/\/3) is decreas- 
ing with the proto-cluster mass since E^, is constant but 

rV^. Thus the more massive clusters can cool more 
efficiently. On the other hand their heating is also more 
intense because it is proportional to M, which scales as 




Fig. 3. Mean gas temperature for various proto-cluster 
masses. 




Fig. 4. Jeans mass vs proto-cluster masses. 



3.3.1. Mean temperature and Jeans mass 

Figure [3] portrays the mean gas temperature as a function 
of proto-cluster mass for a*.c = 1, 2 and 4. Evidently, the 
temperature increases from about 10 K to almost 20 K 
for the biggest mass considered. Since the gas density is 
decreasing as the proto-cluster mass is increasing, this im- 
plies that the Jeans mass is unavoidably increasing with the 
proto-cluster mass as displayed in Fig. 21 Indeed, it shows 
that as expected the Jeans mass within the proto-clusters 
increases from about 5 to ~25 Mq as evolves from 10^ 
to 10^ Mq. This could suggest that, indeed, the IMF could 
vary within clusters. However, as we will see below this is 
not the case. The primary reason is that the peak of the 
IMF also depends i) on the Mach number and ii) on the 
equation of state, i.e. the temperature dependence on the 
density. 



3.3.2. Temperature distribution within proto-clusters 

As discussed in Hennebelle & Chabrier (2009, HC2009), 
the equation of state has a significant influence on the core 
mass function. Knowing the mean temperature only, does 
not appear to be sufficient, therefore. Instead one needs to 
know the complete equation of state, that is to say, how 
the temperature varies with the density within the proto- 
clusters. For that purpose, we compute the temperature at 
various densities, assuming that the turbulent heating Fturb 
and the gradient per km s~^ used for the molecular cooling 
correspond to the mean conditions within proto-clusters of 
mass M^. This means that we are assuming that both the 
heating and the cooling, only depend on the large scale 
conditions and do not vary locally. 



Hennebelle: Star formation within clusters 



9 




log(n} (cm } 

Fig. 5. Temperature vs density for various proto-cluster 
masses. 



where M is the Mach number, 6 ~ 0.5 — 1 (Kritsuk et al. 
2007, Audit & Hennebelle 2010, Federrath et al. 2010) and 
po is the mean density. 

The self-gravitating fluctuations are determined by 
identifying the structures of mass M in the cloud's ran- 
dom field of density fluctuations. These are gravitationally 
unstable at scale i?, according to the virial theorem. This 
condition defines a scale- dependent (log)-density threshold, 
\i\{pc{R) / po)-, or equivalently, a scale-dependent Jeans 



mass 



Mi 



Mi 



2/3 



G 



-R 



iG 



ipcy 



2»7 



R 



(31) 



where Cs is the sound speed, G the gravitational constant, 
oj a constant on order of unity while Vo and 77 determine 
the rms velocity. 



rins/ 



y2 
^0 



i- 

Vlpc 



2>) 



(32) 



Figure [5] shows the gas temperature as a function of 
density for various proto-cluster masses. For densities be- 
tween 10'^ and ~ 10^ cm""^, the temperature decreases from 
about 20 to ~ 8K, leading to an effective adiabatic expo- 
nent of about ~ 0.85 — 0.9. At higher densities, the tem- 
perature increases slightly due to the infiuence of the dust. 
The minimum temperature is reached for n ~ 5 — 8 x 10^ 
cm~^. This compares reasonably well with the tempera- 
ture distributions displayed in Fig. 2 of Larson (1985). The 
proto-cluster mass does not have a drastic infiuence on the 
temperature distribution with variations of a few degrees 
only. 



3.4. Mass distribution: Hennebelle-Chabrier theory 

In this section, we briefiy describe the ideas and the for- 
malism of the Hennebelle-Chabrier theory of star formation 
(HC2008, HC2009, HC2011) that is used in this manuscript 
to infer the IMF within clusters. In this theory the prestel- 
lar cores that eventually lead to protostars and then stars 
correspond to the turbulent density fiuctuations that arise 
in a supersonic medium, which become sufficiently dense to 
be self-gravitating. These fluctuations are counted by using 
the formalism developed in cosmology by Press & Schechter 
(1974) although with a formulation close to the approach of 
Jedamzik (1995). Recently, Hopkins (2012) formulates the 
problem using the excursion set theory (e.g. Bond et al. 
1991) formalism. He also extends the calculation by con- 
sidering the whole galactic disk. His results regarding the 
small scale self-gravitating fiuctuations (the last crossing 
barrier) are undistinguishable from the HC2008 result at 
small masses and only slightly different at high mass. 

According to various simulations of hydrodynamic or 
MHD supersonic turbulence, the density PDF is well rep- 
resented in both cases by a lognormal form. 



ViS) 

s 



1 



: exp 



2al 



-al/2 , al =ln(l + 6^X^ 



,J = ln(p/po), (30) 



Because a fluctuation of scale R is replenished within a typ- 
ical crossing time tr, and is therefore replenished a number 
of time equal to T^j:/{(t)tTB.,ff), where Ti^.jf = T]i/(j)t is the 
freefall time at scale R (see Appendix of HC2011), i.e. at 
density pn Mr/ {Att/3R^), HC2011 includes this con- 
dition into the formalism originally developed in HC2008. 
This yields for the number-density mass spectrum of grav- 
itationally bound structures, Af{M) = d{N/V)/dM 



M{Mr) 



MndMn [ dR^ Wj,'' ^ 



(33) 



which is Eq. (6) of HC2011 and except for the time ra- 
tio is similar to Eq. (33) of HC2008. Note that we 
have dropped the second term which appears in Eq. (33) 
and which entails the derivative of the density PDF. For 
the strongly self-gravitating regimes that we are consider- 
ing here, this approximation is well satisfied. 

Equation ([55]) gives the mass spectrum that we are com- 
puting in the following. It depends on Mn and dMn/dR, 
which can be obtained from Eq. (1311) . To proceed, it is more 
convenient to normalize the expressions. After normalisa- 
tion, Eq. ((3T|) becomes 



(34) 



= M/MO - R i^fip) + MlR^'^) , 
ftp) = f(Mk) ^ 

where C° = Cs{po), M^, Aj and are given by 



M°j = aj 



10^ cm 
aj_y/' (CO) 



2.33' 



(35) 



T 




1/2 /_/f_\-l/2 



10 K 



2.33' 



(36) 
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(0.8-1.0) 



O.lpc 



0.2 km s" 



(37) 



aj and Cm being dimensionless geometrical factors on the 
order of unity. Taking for example the standard definition of 
the Jeans mass, as the mass enclosed in a sphere of diameter 
equal to the Jeans length, we get aj = while Cm = 

47r/3 (e.g. HC2009). 

With Eqs. dSni) and dM]), we obtain 



AA(M) 




(38) 



An important difference to the work of HC2009 is that 
the equation of state shown in Fig. [5] is not polytropic, 
but fully general. To obtain dAIf^/dR, we must differentiate 
Eq. dMl, which leads to (see Appendix of HC2009) 



dM% /(||) - 3^/'(#) + (2ry + l)MlR^^ 



dR 



(39) 



With this last equation, all quantities appearing in Eq. 
are known and the mass spectrum of the condensations can 
be computed. Note that Mfi must be computed numerically 
from Eq. ([M)) . 

4. Results: mass distribution of self-gravitating 
condensations in clusters 

In this section we calculate the mass spectrum of the self- 
gravitating fluctuations and discuss its dependence on the 
various parameters. 

4.1. Preliminary considerations 

Before presenting the complete distribution, we start by 
discussing the position dependence of the distribution peak. 
For that purpose we make various simplifying assumptions. 

4.1.1. Peak position: dependence on the Jeans mass and 
Mach number 

Unlike what is often assumed in the literature, it is un- 
likely that the peak of the CMF/IMF solely depends on 
the Jeans mass. In particular, it likely depends on the Mach 
number because compressible turbulence creates high den- 
sity regions where the Jeans mass is smaller. Indeed, in any 
turbulent medium, there is a distribution of Jeans masses 
rather than a single well-defined value. The peak position 
has been calculated by HC2008 who show (their Eq. 46 ) 
that Mpeak = My{l + b^M^f/* (note that in HC2008 b 
stands for b^ as different notations were used). However, 
as discussed above, time-dependence is not considered in 
IIC2008 while it is taken into account in Eq. (p8| through 
the term I/th. To calculate the peak position in this case, 
we proceed as in HC2008, that is to say we neglect the 



turbulent support (i.e. we set =0), which has little in- 
fiuence on the peak position. We also assume strict isother- 
mality within the cluster, that is to say, we set / = 1. As dis- 
cussed in IIC2009, the equation of state indeed has an influ- 
ence on the peak position. However, as portrayed in Fig. [5l 
the deviation from the isothermal case is not very impor- 
tant with an effective adiabatic index of 0.85-0.9. Moreover, 
our goal here is more to discuss the dependence qualita- 
tively rather than getting an accurate estimate, which is 
calulated later in the manuscript. 

Assuming A^* = and / = 1, it is easy to take the 
derivative of Eq. psp and to show that the maximum of 
M(M) is reached for M = 1/(1-1- 6^A^^), which implies 
that the distribution peak is reached at 



peak - ^ ^ ^/^2 ■ (40) 

The difference to the expression presented in HC2008 comes 
from the fact that more small structures form when time 
dependence is taken into account since the small scale fluc- 
tuations are rejuvenated many times while the larger ones 
are still evolving. 

Note that as already stressed, the peak position does 
not depend on the Jeans mass only, but also on the Mach 
number. For a typical b^ of the order of ~0.5 (see Federrath 
et al. 2010) and a Mach number of about 5, we derive that 
the peak position is typically shifted by a factor of about 
10 with respect to the mean Jeans mass. 



4.1.2. Dependence of the peak position on the cluster mass 

To estimate the peak position, we must therefore estimate 
the Jeans mass and the Mach number as a function of 
the cluster parameters. Since the Jeans mass is equal to 
(7r5/76)C3/Vp^;G, we obtain that 

Mpeak ^ 6 62 Gl/2(p,) 1/2^2- (^1) 

To estimate the sound speed we must compute the mean 
temperature within the proto-cluster. For the sake of sim- 
plicity, we consider in this analytical estimate the turbu- 
lent heating Tturb and the molecular cooling Kmoi given by 
Eqs. (PT|) and (1^ . respectively. Writing Tturb = ^moi/n, 
we obtain the temperature and thus the sound speed 



C = — 
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On the other hand, with Eqs. ([T 
1 
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Ipc^ 
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Fig. 6. Expected position of the CMF peak as a function 
of the proto-cluster mass for 3 values of the ratio of parent 
clump over proto-cluster mass ratio, a*,c. 



As already mentioned, typical values of d are 2-2.5 while 
c is typically low and close to zero. Ignoring the depen- 
dence on n'^, the peak position consequently has a weak 
dependence on the cluster radius, which typically depends 
on R as RI^^~^. This is the result of partial compensation 
of the sound speed, the density and the velocity dispersion. 
Thus we can conclude that in the regime where turbulent 
heating dominates over the cosmic ray heating, the initial 
mass function is expected to weakly depend on the cluster 
masses. According to our estimate, this is the case for large 
mass clusters, i.e. with masses above ~ IO^Mq. Note that it 
should be kept in mind that in this analysis, isothermality 
is assumed, which is also a simplification. 

Before solving for the whole mass spectrum obtained 
with a complete equation of state, it is worth comput- 
ing the dependence of Mpcak — + h'^M'^) on the 
proto-cluster mass and on Qq cx a^'^^, the parent clump 
over proto-cluster mass ratio using the complete ther- 
mal balance stated by Eq. rather than the simpli- 
fied approach used above. Figure ([6]) portrays the value 
of Mpeak = Mj/{1 + b'^M^) as a function of the proto- 
cluster mass for three values of a*.c and b — 0.6. This last 
value corresponds to a mixture of solenoidal and compres- 
sive modes (Federrath et al. 2008). The trends are sim- 
ilar to what has been inferred from Eq. (pS)) . that is to 
say, the peak position weakly depends on the proto-cluster 
mass and more strongly on the accretion rate controlled by 
a*_c- However, it is remarkable that for proto-clusters whose 
masses are between 10^ and 10^ Mq, and within the ranges 
Q!*,c = [1,4], which encompasses most of the data points 
reported in Fig. [21 the peak position varies by a factor of 
only 3 and only 2 for between 10"^ and 10^ M0. 



4.2. Fiducial parameters 

We now present the complete mass spectrum obtained by it- 
eratively solving Eq. and computing Eqs. and ([55)1 
using the temperature distribution obtained in section l5.3.2l 
and portrayed in Fig. [5] Because the latter is not analytical, 
it is difficult to take its derivative as required by Eq. (15^ . 
Thus for this purpose we first obtain a fit of the tempera- 
ture using a high order polynomial. 

Upper panel of Fig. [7] portrays the resulting mass spec- 
tra for various cluster masses as labeled in the figure. The 
dots represent the Chabrier IMF (Chabrier 2003) shifted by 
a factor of about 2. Indeed, since we compute the mass func- 
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Fig. 7. Mass spectra for various cluster masses and the fidu- 
cial value a*_c = 2. The dots represents the Chabrier IMF. 
The two lines represents the distribution and M~^ '^. 
The upper panel displays the mass distribution obtained 
with Eq. (pS)) while the lower panel shows this mass distri- 
bution convolved by a Gaussian of dispersion a = M /2. 



tion of self-gravitating condensations, this shift is necessary 
to account for the observed shift between the CMF and the 
IMF (e.g. Alves et al. 2007, Andre et al. 2010) and for the 
theoretical estimate for the core-to-star efficiency (Matzner 
& McKee 2000, Ciardi & Hennebelle 2010). Because we do 
not consider any cluster mass distribution and efficacity 
problem here, we set the distribution maximum to 1. This 
allows us to compare the shape of the various distributions 
more easily. 

Various points are worth discussing. First of all, the dis- 
tributions all agree with the Salpeter and Chabrier IMF at 
high masses. As discussed in HC2008, this is because the 
high mass part is essentially due to the turbulent support, 
which does not change much in this regime. Moreover, the 
high mass part of the distribution tends to slowly depend 
on the parameters that control the turbulence, namely M 
and M^. Second, the peak position varies by about a factor 
2 as the proto-cluster mass changes from ~ IO^Mq to about 
lO'^M© while the distributions corresponding to M, = 10^ 
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and 10'' Mq are nearly indistinguishable. This weak varia- 
tion is, as discussed in the previous section, a consequence of 
various quantities compensating each other. Third, the po- 
sition of the peak itself is close to the peak of the Chabrier 
IMF shifted by a factor of 2. 

4.3. Link between the CMF and the IMF 

An important question regarding the mass distribution is 
the link between the IMF and the CMF, which is not ex- 
pected to be as simple as a unique efficiency (e.g. Alves et 
al. 2007, Goodwin et al. 2008) on the order of 2-3. It is 
instead expected that the two distributions are correlated 
with some dispersion because the initial mass and veloc- 
ity distributions within the cores vary from one core to 
another, therefore leading to different evolutions. In other 
words, the virial theorem is not representing the problem in 
its full complexity and the initial conditions and the bound- 
ary conditions (which enter through the surface terms, see 
e.g. Dib et al. 2007) influence the mass of the objects that 
form within the collapsing cores. Indeed, Smith et al. (2008) 
have shown that in hydrodynamical simulations, the corre- 
lation between the mass of the sink particles and the cores 
in which they are embedded is very good during the first 
freefall times and becomes less tight after ~ 5 freefall times. 
This is because initially the sink particles are accreting the 
mass of the parent core. As time evolves, the material that 
falls onto the sinks comes for further away and is less and 
less correlated to the initial mass reservoir. More quantita- 
tively, Chabrier & Hennebelle (2010) have shown that in the 
simulations of Smith et al. (2008), the correlation between 
the CMF and the sink particles mass function, which is 
likely representing the IMF, can reasonably be reproduced 
at intermediate time (say 3-5 freefall time) by a Gaussian 
distribution with a width, a, on the order of M/3—M/2. To 
take this into account, we have convolved the mass distribu- 
tion portrayed in the upper panel of Fig. [7] with a Gaussian 
distribution of width aconv = M/2. 



— 

dM conv 



rS(*^')exp(-^^) 



dM' 



Jo cxp 



-{M-M') 



dM' 



(45) 



Note that in principle, since the efficiency of the accretion 
is lower than one, the distribution should be shifted toward 
smaller masses (e.g. Alves et al. 2007, Andre et al. 2010). 
In practice, to facilitate the comparison with dN/dM, we 
have not shifted it here. 

As expected, the distribution dN/dMconv is slightly 
broader than the distribution dN/dM in particular at low 
masses whereas, the high mass part and the peak position 
are almost unchanged. This effect could therefore in partic- 
ular contribute to the formation of low mass brown dwarfs. 
To clarify, because some cores are marginally bound, for ex- 
ample because their velocity field is initially globally diverg- 
ing, few objects significantly less massive than the parent 
core mass form, most of the envelope being then dispersed 
in the surrounding ISM. 

To summarize, for an accretion rate that allows one 
to reproduce the mass-size relation of clusters, we self- 
consistently predict a distribution of self-gravitating struc- 
tures that i) is very close to the field IMF inferred by 
Chabrier (2003), ii) is almost independent on the cluster 
mass for gas masses larger than lO'^Af©. We stress that if 



more variability is expected for low mass clusters, it is very 
difficult to infer any reliable statistics and thus no data are 
available in this regime. 

4.4. Dependence on the accretion rate 

Because the accretion rate onto the proto-clusters of a given 
mass is likely varying, as suggested by Fig. [2l we explore 
the effect of this parameter on the whole mass spectrum. As 
anticipated in Fig. El it has some infiuence on the peak po- 
sition. As discussed in section 2.1, we use the parent clump 
mass over proto-cluster mass ratio to quantify the accre- 
tion rate following Eqs. dS])-® and we select the two values 
a*_c — 1 and 4, which as shown by Fig. [2] encompass almost 
all points of the observational distribution. 

The mass distributions obtained for a^^c — 1 and 4 
are displayed in upper and lower panels of Fig. |51 respec- 
tively, which shows that the trends inferred in Fig. [5] are 
well confirmed. The mass distribution is shifted toward 
larger masses for a lower accretion rate (a*_c = 1) and to- 
ward smaller masses for a higher accretion rate (a*,c = 4). 
Apart from this, the general behaviour is almost identical 
to the fiducial case a*.c = 2, that is to say, the distribu- 
tions weakly depend on the proto-cluster gas masses for 

~ 10'^ — 10^ M0 and the high mass part is always close 
to the Salpeter IMF. 

Importantly enough, we note that within the range 
1 < a*_c < 4 the variation of the mass spectrum with the 
accretion rate remains limited to a factor of about 2 almost 
everywhere, especially for the proto-clusters, whose mass of 
gas is between 10"^ and 10^ Mq. 

4.5. Dependence on cosmic rays and radiation 

Another source of possible variability are the heating 
sources, namely cosmic rays and interstellar radiation field, 
which heat the gas and the dust respectively. In this section, 
we investigate the dependence of the CMF on these param- 
eters. Indeed, both are likely varying as one approaches a 
supernova remnant or a massive star for example. 

4.5.1. Dependence on cosmic rays 

Figure [5] shows the mass distribution for two values of cos- 
mic rays heating, namely half (upper panel) and three times 
the fiducial value stated by Eq. (|23)) . As can be seen from 
a comparison with the upper panel of Fig. [71 cosmic ray 
heating has an impact on the mass distribution that is only 
moderate. For most of the distribution, the variations are 
typically less than a factor two. 

4.5.2. Dependence on radiation 

When increasing the interstellar radiative fields, the var- 
ious components must be distinguished since they have a 
different physical origin. For simplicity we multiply all com- 
ponents by the same factor except for the cmb component, 
which is kept constant. Figure [10] shows the mass distri- 
bution for two values of interstellar radiation field, namely 
half (upper panel) and twice the fiducial value. The depar- 
ture from the mass distribution displayed in Fig. [7] is even 
smaller than for the cosmic rays. This is because the ra- 
diation field is modifying the dust temperature, which has 
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Fig. 8. Mass spectra for various cluster masses, upper 
panel a*,c = 1, lower panel a*_c = 4. As in Fig. [7] the 
dots represents the Chabrier IMF. The two lines represents 
the distribution M"^ and M"^-^. 

only a weak impact on the mass distribution because, as 
discussed above, the effective 7 at low and high density are 
about 0.85-9 and 1, respectively (see Fig. [5]). 

4.6. Discussion 

One of the central questions regarding the IMF is that of 
its apparent universality, that is to say, its apparent lack 
of variations in the different determinations that have been 
made so far (e.g. Moraux et al. 2007, Bastian et al. 2010). 
The various estimates in open clusters and young clusters 
lead to some variability on the IMF parameters, in par- 
ticular on the characteristic mass, or peak position, but 
as pointed out by de Marci et al. (2010) dynamical evolu- 
tion is a plausible explanation to account for the strongest 
variations. There is a possible variation independent of the 
cluster dynamical evolution, however, as can be seen for ex- 
ample in Fig. 3 (left panel) of Bastian et al. (2010) and in 
Fig. 2 of de Marci et al. (2010), which shows that at a given 
dynamical age, the value of rric typically varies over a factor 
of about 2. It is worth stressing that a possible deficiency of 
brown dwarfs is observed in Taurus (e.g. Guieu et al. 2006, 
Luhman et al. 2009, Bastian et al. 2010) while a possible 
excess (factor of about 2) may have been observed in a Ori 




log(M) (M,) 



Fig. 9. Same as Fig.[7]for the fiducial value a,,c = 2 but for 
two values of cosmic rays. The upper panel shows cosmic 
ray heating twice lower than the fiducial value while the 
lower panel shows cosmic ray heating three times higher. 

and Upper Sco (e.g. Caballero et al. 2007). These variations 
need to be confirmed, however. It is clearly difficult, at this 
stage, to give a definite conclusion. 

On the other hand, the results obtained in this paper 
can be broadly summarized as follows: i) the high mass 
part of the IMF is robust to all variations explored in this 
paper (higher accretion rate could lead to some flattening 
and shift the peak toward smaller masses), ii) the low mass 
part varies with the mass of the cluster and the accretion 
rate by a factor of about 2 in the range we explored. The 
variation with the cluster mass in the range 500-10* Mq 
is very limited, typically less than a factor 2 (see the four 
more massive cases displayed in Fig. [T]) and thus account 
well for the lack of variability since as displayed in table 1 
of de Marci et al. (2010), most of the clusters that we stud- 
ied have a mass of stars that is comparable to or greater 
than, 10^ Mq. This in turn corresponds to a mass of gas 
about five times larger, assuming an efficiency of 20%. The 
variation of the mass function with the accretion rate, i.e. 
with Q!*^c, is more significant but remains compatible with 
the range of values that we inferred from the comparison 
with the data from Lada & Lada (2003), that is to say a 
factor of about 2 for the accretion rate (at fixed mass) and 
a factor of about 2-3 on the low mass part of the distribu- 
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Fig. 10. Same as Fig. [7] for the fiducial value a*^c = 2 but 
for two values of the radiation field. The upper panel shows 
an interstellar radiation field equal to half the canonical 
values while the lower panel shows a radiation field twice 
as high. 

tion (with much stronger variations for the very low mass 
objects). 

5. Conclusion 

We have developed an analytical model that successfully re- 
produces the mass-size relation of embedded clusters. It is 
based on the continuous accretion of gas that drives the tur- 
bulence which in turn resists the gravitational contraction 
and determines the cluster radius. Comparing the mass-size 
relation with the available data of embedded clusters, we 
showed that they can be well fitted by employing an accre- 
tion rate that is entirely reasonable and close to the rate 
that can be inferred using Larson relations. Moreover, the 
variation of this parameter needed to reproduce the dis- 
persion of observational data is moderate being equal to 
a factor of about 2. Eventually, the turbulent energy dis- 
sipates and heats the gas. Performing a thermal balance, 
we calculated the temperature distribution within the clus- 
ter and applied a time-dependent version of the HC2008 
theory to obtain the mass spectrum of self-gravitating con- 
densations. The peak position and more generally the whole 
mass spectrum is inferred from a large domain of proto- 



cluster masses. We found that for gas masses between 500 
to 10"* Mq corresponding to a mass of stars roughly 5 
times smaller, the mass spectrum of self-gravitating con- 
densations does not vary significantly. The mass spectrum 
varies more significantly with the accretion rate, but given 
the range of accretion rate deduced from the observational 
comparison, they remain compatible with the available de- 
termination of the IMF in the Galaxy. 
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Appendix A: Virial theorem for an accreting system 

The virial theorem is derived in many textbooks (e.g. Shu 
1992, Kulsrud 2005, Lequeux 2005) and does not need to 
be presented. However, it is most of the time assumed that 
the system has a constant mass, i.e. is not accreting. Here 
we derive its expression for a system accreting at a rate M 
(see also Goldbaum ct al. 2011). That is to say, at radius 
R{t), gas at density Pinf is falling into the cluster with a 
velocity Vinf. 

We start as usual by multiplying the momentum con- 
servation equation by the vector position r and integrating 
over the volume. The right-hand side of the momentum con- 
servation equation is equal to the sum of the forces, namely 
thermal pressure and gravity (and Lorentz force when the 
magnetic field is considered). It can be rewritten in the un- 
magnetized case as 2Etherm-^PthermV+Eg, where Etherm 
is the thermal energy, Ptherm is the thermal pressure, V is 
the volume and Eg the gravitational energy. The left-hand 
side is given by / pdtViTidV . Integrating this expression by 
part, we have 



j pdiVindV = J pdt{viri)dV - J pvfdV, 



j pdt{viri)dV + J p-v.grad{viri)dV{A.l) 

dt{ j pviridV) + j diy{pViriv)dV 
- 1 pv^dV, 

dt{ J pdt{r'l)/2dV) + j pVin^rdS 
- j pvldV, 



where we used the continuity equation dtp + div(pv) = 
to cancel some of the terms. 

Let us consider now the first term which appear in the 
right-hand side. We can write 



/ 



pViudV = 



j pdt{r^)/2dV = J pdt{r1)/2dV 
+ j /9V.grad(rf/2)ciy, 
= J dt{prf/2)dV + J div{prfv)/2dV, (A.2) 
= dtij pry2dV) + j pr^/2v.dS. 
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Thus Eqs. (|X2|) and (|X3]) lead to 
j pdmrdV = df{ J prf /2dV) + dt{ J prf /2^.dS) 
pviVivdS — / pv'fdV, 



l/2dt.I-dt{pv,^fR'S) (A.3) 



+Pinfv1nfRS - 2E, 



= l/2dll - \dt{M)R^ ~ \M^dMR'' 

where / = / prfdV, M = Pi„fVinfS and Pram = Pinf^l^f 
Note that the minus sign in the second and third terms of 
the right-hand side, are due to the fact that as the gas is 
infalling, Vinf < and thus v.dS < 0. 

The complete expression for the virial theorem is thus 

l/2d^2l = ^dtiM)R^ + ^M^dMR^ - KPram + Ptherm)V 
+2E, + 2Etherm+Eg. (A.4) 

Assuming stationarity and neglecting thermal support 
and thermal pressure, we derive 



-M^OmR^ - iPrarnV + 2E, + Eg = 0, (A.5) 



1 

2 

which is the equation used in this paper 



Appendix B: Gravitational energy released by 
accretion 

We calculate here the gravitational energy that is released 
when a particle of gas is accreted onto the proto-cluster. 
For that purpose, we assume the following. First, the ac- 
creted particles that enter the proto-cluster, will on average 
be distributed uniformly within the cluster, second, we as- 
sume that the gravitational energy inside the cluster is pro- 
duced by a uniform distribution of matter, i.e. we neglect 
the influence of any local density fluctuations. 

Under these assumptions the gravitational force inside 
the proto-cluster is — G47r/3/9*r and the gravitational po- 
tential is thus G27r/3(r2 - 3Rl). Since the gravitational 
potential outside the cluster is —GM^/r, it can easily 
be checked that the gravitational potential is continuous 
throughout the cluster boundary, i.e. at r = i?*. Therefore 
we get that the mean specific gravitational energy, eg, of a 
particle that falls into the proto-cluster is 

— ^ — (3i?* - r )f{r)dr = -5-^, (B-1) 

where f{r)dr — 3r^dr/Rl is the fraction of particles located 
between r and r + dr. 

Thus the total accretion energy rate is 

5 K* 



